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Abstract 

Color-octet contributions to the associate J/ip+j are found to be negligible, 
compared to the ordinary color-singlet contribution. Within the color-singlet 
model the J/ip + 7 production in the leading order is possible only through 
gluon-gluon fusion process. Therefore, the associate J/'0+7 production remains 
to be useful as a clean channel to probe the unpolarized and polarized gluon 
distribution inside proton and to study heavy quarkonia production mechanism. 
We discuss in detail the associate production of J ftp + 7 at pp (or pp) and ep 
colliders. By requiring the J ftp to decay into an e + e~ or u + u~ pair, we end up 
with an exceptionally clean final state. This process can therefore serve as a 
very clean probe of the gluon densities inside the proton as well as the photon. 
Numerical results are presented for the TEVATRON pp and HERA ep colliders. 
This same mechanism can be used to probe the polarized gluon content of the 
proton in polarized p + p(p) collisions. We study in detail J ftp + 7 production 
at both polarized fixed target and polarized collider energies for RHIC. 



1 Introduction 

The production of heavy QQ bound states (quarkonia) offers a good testing ground 
for perturbative QCD, since here one combines relatively large cross sections with 
rather clean final states. The large cross sections are, of course, due to the strong 
(colour) interactions of the quarks Q; clean signals emerge when the s-wave vector 
states (J/ip or T) decay into a pair of charged leptons (e + e _ or /i + fi~). Unfortunately 
the theoretical analysis of the inclusive production of quarkonia in hadronic collisions 
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is complicated, due to the large number of contributing processes. For instance, the 
following processes contribute 0] to hadro-production of J/ip: 



9 + 9 






9 + 9 ' 


Xi(~ 


-> J/V> + 7) + g\ 


g + q - 


Xi(- 


-> J/^ + 7 ) + g; 


q + q - 


x*(- 


-> J/^ + 7) + <?; 


g + g,q + q 


-&(- 


> J/^ + X) + 6. 



For the photo- (or lepto-) production of J/if), one has to consider the processes 

i + g ^J/*P + g; 

j + g ^b(^ J/^j + X) + b. (1.2) 

At the relatively low energies of fixed target experiments, only the first reaction in 
( p.. 2D leads to a sizeable cross section. Indeed, recently the NMC collaboration has 
shown |J that this reaction can be used to determine the large- x behaviour of the 
gluon distribution inside the proton. However, at the much higher energies that can 
be achieved at the HERA collider the situation is considerably more complicated; 
here one does not only have to include || J/ip production from b decays, but also all 
processes of (|1.1|), since at these energies the quark and gluon content of (quasi- 
)real photons can no longer be ignored. Indeed, J/ip production at HERA has been 
suggested || [Tj as a probe of the gluon content of the photon. At the same time, one 
hopes to constrain || the small- x behaviour of the gluon content of the proton using 
reactions (|1.2j ). In order to achieve both goals, one will clearly have to discriminate 
between the different mechanisms of J/ ip production, and various methods to do this 
have been suggested [0. 

Another staple of perturbative QCD is the production of direct photons in 
hadron-hadron collisions |§ ; more recently, these calculations have also been extended 
to ep colliders pi. In leading order, direct photons can be produced by gq scattering 



or qq annihilation, so the analysis of such events at high energy pp or ep colliders is 
again quite complicated. The situation is simplified when one requires |TTJ a heavy c 



or b quark to be produced together with the photon, but these reactions are not well 
suited to measure gluon structure functions, since the cross sections also depend on 
the poorly understood heavy quark distribution functions. This difficulty is avoided 
if we move the second heavy quark from the initial to the final state; if we in addition 
require the two heavy quarks to form a bound state, we can expect a very clean final 
state, as discussed above. In this paper we therefore study the associate production 
of a QQ bound state and a hard, isolated photon. We focus on the production of J / tp 
states, which offer the largest rates, but the generalization to other s-wave J = 1 cc 
or bb states is straightforward. In leading order the J /ip + 7 final state can only be 



produced in gg fusion [12 



g + g ^J/^ + 7; (1-3) 

g + g ^Xc^J/^ + i- (1.4) 
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The reaction (O ) can only produce photons with small transverse momentum p 7 



J T 1 

since we are interested in the production of isolated photons with high p T , we only 
need to consider reaction ( |1.3|) . (Note that g + g — > Xc + 1 is not possible.) The cross 
section is then given by 

da(A + B^J/ijj + 1 + X) = J dx 1 dx 2 f g \ A (x 1 )f glB (x 2 )da(g + g -f J/ijj + 7), (1.5) 

where A and 5 can be a hadron, photon, or electron. We will use the so-called colour 
singlet model to estimate the hard subprocess cross section a. In this model the J/ip 
is treated as a nonrelativistic cc bound state. One then has [P"3]| 

d8j£+g-^J/f+j} _ 16Traa 2 m^,\R(0)\ 2 
dt ~ 27s 2 



t _|_ u 

(u-ml) (S-ml) (s-rnl) (t-mj) 



(1.6) 

here, s, t and u are the Mandelstam variables of the parton-parton collision, = 
3.1 GeV is the mass of the J/tp meson, and the cc wave function at the origin |i?(0)| 2 
can be determined from the leptonic decay width of J/ip, in leading order expression: 

T(M -> e+e") = — ^r|i?(0)| 2 = 4.72 keV => \R(0)\ 2 = 0.48 GeV 3 . (1.7) 

2 Color-Singlet and Color-Octet Contribution in 
J/^ + 7 Production 

[For the details of this Section, see Ref. fi~4]| .l Up to very recently the conventional 



color-singlet model [J3| had been used as only possible mechanism to describe the 
production and decay of the heavy quarkonium such as J/ip and T. An artificial 
.fT-factor (~ 2-5) was first introduced to compensate the gap between the theoretical 
prediction and experimental data. The uncertainties related to this large i^-factor 
are the following: We don't know precisely the correct normalization of the bound 
state wave function, the possible next-to-leading order contributions, the mass of the 
heavy quark inside the bound state, etc. However, even with this large i^-factor some 
experimental data are found to be difficult to describe. For the direct J / ip production 
with large p T in pp collisions, the dominant mechanism has been found to be through 
final parton fragmentation [ITS] . As applications of this fragmentation mechanism, 
various studies of prompt charmonium production at the Tevatron collider |16|, |17j 
have been carried out, and the CDF data on the prompt J/ip production ]T9 



qualitatively meet these theoretical predictions. Nevertheless, the if)' production rate 
at CDF is about 30 times larger than the theoretical predictions, which is the so- 
called ip' anomaly, even after considering the fragmentation mechanisms of g — > ip' 
and c — > ip'. As a scenario to resolve this ip' anomaly, the color-octet production 
mechanism was proposed ||20|| . By using this idea, heavy quarkonium (charmonium) 
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hadropro duct ion through the color-octet (cc)s pair in various partial wave states 
( 2S+1 Lj) has been studied in addition to the color-octet gluon fragmentation approach 
pl|, |22|. We note that the inclusive T production at the Tevatron also shows an excess 
of the data over theoretical estimates based on perturbative QCD and the color-singlet 
In this case, the p T of the T is not so high, so that the gluon fragmentation 



model 23 



picture may not be a good approximation any more. 

Since it has been proposed that the color-octet mechanism can resolve the if) 
anomaly at the Tevatron collider, it is quite important to test this mechanism at 
other high energy heavy quarkonium production processes. Up to now, the following 
processes have been theoretically considered: inclusive J /if) production at the Teva- 
tron collider and at fixed target experiments [21 , ^2| [24fl , spin alignment of leptons 
to the decayed J /if) [f^|, the polar angle distribution of the J /if) in e + e~ — > J /if) + X 
26fl , inclusive J /if) productions in B meson decays [27, |28| and in Z° decays at LEP 
, J /if) photoproduction ]28|, |33[] , and color-octet J /if) production in 



T decays p5[ . For more details on recent progress, see Ref. 
the theoretical developments on the quarkonium production 



which summarizes 



Recently, the helicity 
decomposition method in NRQCD factorization formalism was developed by Braaten 
and Chen [[37]]. With this method, polarized J /if) production in B decay was consid- 

H 



ered in Ref. 
of different 3 



381. In Ref. [37 



it was pointed out that there are interference terms 
Pj contributions in polarized J /if) production. 
In this Section, we investigate the color-octet mechanism for associate J /if) + 7 
production in hadronic collisions. The associate J /if) + 7 production has been first 
proposed as a clean channel to probe the gluon distribution inside the proton or 
photon fl^| , and then to study the heavy quarkonia production mechanism [[5], f|l 
as well as to investigate the proton's spin structure [J42|| . If the color-octet mechanism 



gives a significant contribution to this process, the merit of this process to probe the 
gluon distribution etc. would be decreased or one should find suitable cuts to get rid 
of this new contribution. Associate J/if> + 7 production has also been studied as a 
significant QCD background to the decay of heavier P-wave charmonia (x c j{P)) in 

53 1 within the color-singlet model 
(c or g) + 7 + X 



fixed target experiments 



In Ref. [g§,the 
J /if) + 7 + X at Tevatron 



fragmentation contribution from p + p 
energies ( v /s=1.8 TeV) is studied, and these fragmentation channels are found to be 
significantly suppressed compared to the conventional color-singlet gluon-gluon fusion 
process. 



2.1 General Discussions 

Now we consider the associate production of a J /if) and a photon at hadronic colliders 
and fixed target experiments, including the new contribution from the color-octet 
mechanism in addition to the usual color-singlet contribution. The only possible 
subprocess in the framework of the color-singlet model in the leading order is through 
gluon-gluon fusion, 

g + g - 7 + (cc)( 3 ^ 1} )H JM- (See Fig. 1) (2.8) 
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This gluon-gluon fusion process for the associate J /if) + 7 production is again possible 
within the color-octet mechanism. (See Fig. 1). However, in the leading order color- 
octet contributions there also exist many more subprocesses, which are 

9 + 9 - 7+(cc)C4 8) )H J/1>), (2-9) 

9 + 9 - 7 + (cc)( 3 Pf)H<m (2-10) 

q + q - 7 + (cc)( 1 4 8) )H W), (2.11) 

g + g - 7 + (cc)( 3 Pf)H W), (2.12) 

q + q - 7 + (cc)( 3 ^ 8) )H W (2-13) 

These subprocesses are possible through the effective vertices of 

7 + S->(cc)( 1 ^ 8) or 3 pf), (2.14) 
which is shown in Fig. 2, and 

q + q^(cc)Cs[ 8) ), (2.15) 



which is shown in Fig. 3. We note that the ( 3 5'| 8 ^) channel in Eq. ( |2.14|) is vanishing 



even though the gluon is not on its mass shell. The Feynman diagram for the gluon 
initiated subprocesses, Eqs. ( |2.9| , |2.10|) , is shown in Fig. 4 (a). For the case of quark 
initiated subprocesses, there are two kinds of diagrams; The Feynman diagram for 
quark initiated subprocesses Eqs. (p. 11] , |2.12| ), due to the effective vertex Eq. (|2.14j ), 



is shown in Fig. 4 (b). The other diagrams, which are possible through the effective 
vertex, Eq. ( |2.15|) , are shown in Fig. 5. 

We first perform a naive power counting analysis of those various subprocesses. 
Free particle amplitude-squared for the color-singlet gluon-gluon fusion process (|2T8| ) 
is [C^aa^)]. In the NRQCD framework [4~7|, the color-octet (QQ) states can also form 



a physical J/ip state with dynamical gluons inside the quarkonium with wavelengths 
much larger than the characteristic size of the bound state (~1/ (MqVq)). In Coulomb 
gauge, which is a natural gauge for analyzing heavy quarkonium, these dynamical 
gluons enter into the Fock state decomposition of physical state of J/ip as 

|JM = OWKQQX 3 ^)) + 0(v Q )\(QQ)( 3 PP)9) 

+ 0(v 2 Q )\(QQ)( 3 s[ h8) )99) + 0(v 2 Q )\(QQ)CS^)g) + ... . (2.16) 

In general, a state (QQ)( 2S+1 Lj) can make a transition to (<5<3)( 25+1 (£ + 1)j)> or 
more specifically (QQ)( 3 p( 8) ) -> (QQ)( 3 ^ 1} ) through the emission of a soft gluon 
(chromo-electric dipole transition), and it is an order of vq suppressed compared to 
the color-singlet hadronization. For the case of chromo-magnetic dipole transition, 
such as (QQ)i 1 So 8 " 1 ) — ► (QQ)( 3 Si^), it is suppressed by Vq at the amplitude level. If 
we also consider the fact that the P wave state is MqVq order higher than S wave 
state in the amplitude level, one can naively estimate that the transitions 

(gg)( 3 pj 8) ) - J/^j + X and (QQ)C4 8) ) -> + X (2.17) 
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are commonly Vq order suppressed at the amplitude-squared level, compared to the 
transition 

(QQ)( 3 Si 1] )^ J/i). (2.18) 

Since these color-octet processes have the same order 0(aa 2 s ) as that of the color- 
singlet gluon-gluon fusion process in free particle scattering amplitude, the color-octet 
subprocesses are order Vq suppressed compared to the color-singlet gluon-gluon fusion 
subprocess. 

However, such analyses are only naive power counting, and do not guarantee 
that the color-octet contributions are suppressed compared to the color-singlet gluon- 
gluon fusion contribution all over the allowed kinematical region. If we consider 



inclusive J/ip photoproduction via 2 — > 2 subprocesses, as shown in Refs. |28], [32 
color-octet contributions dominate in some kinematical region, even if the naive power 
counting predicts the suppression of color-octet contributions compared to the color- 
singlet one. In this respect, it is worthwhile to investigate in detail how much the 
color-octet mechanism contributes to the J/ip + 7 hadropro duct ion. 

In order to predict numerically the physical production rate, we need to know 
a few nonperturbative parameters characterizing the fragmentation of the color-octet 
objects into the physical color-singlet J/ip. These nonperturbative matrix elements 
for the color-octet operators have not been determined completely yet. After fitting 
the inclusive J/ip production at the Tevatron collider, using the usual color-singlet 
J/ip production, the cascades production from x c (lP), and the new color-octet con- 
tributions, the authors of Ref. p2| have determined 



(0|O|( 3 5i)|0) = (6.6 ±2.1) x 1(T 3 GeV 3 , (2.19) 

msm + m^m = (2 , 2±05)xl(r2GeVSj (220) 



c 



with M c = 1.48 GeV. Although the numerical values of the above two matrix elements, 
(0|^( 3 P )|0) and (0\Ot ( x S )\0) , are not separately known in Eq. ( ggg ), one can still 
extract some useful information from them. Assuming both of the color-octet matrix 
elements in Eq. (12.20 ) are positive definite, then one has[^8|] 



< (0\OtCS )\0) < (6.6 ± 1.5) x 1(T 2 GeV 3 , (2.21) 
< (°\°t^o)\0) < {22 ± 5) x 1Q -2 GeV s_ {2 22) 

c 

These inequalities could provide us with a few predictions on various quantities related 
to inclusive J/ip productions in other processes, and also enable us to test the idea 
of color-octet mechanism in the associate J/ip + 7 production process. 

3 The matrix element (0|Og ( 3 Po)|0) can be negative whereas (0|Og( 1 5o)|0) is always positive 
definite Q . Here we choose these ranges just for simplicity. 
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2.2 Kinematics 

Let us consider the process 

a(Pi)Mft) + KP2) M p 2) - JMP) + 7(*)- (2-23) 

in pp (or alternatively pp) collisions. We can express the momenta of the incident 
hadrons (p, p) and partons (a, b) in the pp CM frame as 

P 1 = 2^(1, +1,0), Pl = |(xi,+ii,0), (2.24) 

P 2 = ^(1, -1,0), p 2 = ^(x 2 , -x 2 , 0), (2.25) 

where the first component is the energy, the second is longitudinal momentum, and 
the third is the transverse component of the particle's momentum. The variables x\ 
and x 2 are the momentum fractions of the partons. The momenta of the outgoing 
particles are given by, 

P = (E^,P?,+P T ) = (M T coshy^,M T smhy^,+P T ), (2.26) 
k =(E\P2,-P T ) = (P T coshy\ P T smhy\-P T ), (2.27) 

where Pt is the common transverse momentum of the outgoing particles, Mt is the 
transverse mass of the outgoing J ftp, and y^ (or y 1 ) is the rapidity of J/ip (or 7). 
In order to get the distributions in the invariant mass Mj/^ +1 {= \fs) and transverse 
momentum P T for the process g+g — > J/^+7 process, after introducing dimensionless 
variables, 

x T = 2P T /v/s, x T = 2M T /Vs and r = M|/s, (2.28) 
we express the differential cross section as, 

d<J = f g /p(xi,Q 2 )f g /p(x2,Q 2 )^dx 1 dx 2 di 

= fg/ P (xi, Q 2 )f g /p(x 2 , Q 2 )^g J y XlXT M J/lp J dxidx T dMj/^ +1 

= f 9 / P (xu Q 2 )Mx2, Q 2 )%J (^f) d Xl dy^dP T , (2.29) 

at \Xiy v fTl 



where the corresponding Jacobians are given by 



j ( xix 2 i \ _ 2x 2 x T Mj/^ +1 j { x x x 2 i \ _ A.x x x 2 P T 

\x x x t Mji^) ~ x T (x 2 e+y^ -xtf-v*) an [x^PtJ ~ 2x 1 -x T e+v i 



(2.30) 



Then the distributions are expressed as 



<h f d Xl dx T J ( X \ X f ) ^ , , (2.31) 



dMj/^ +1 J ' \xix T Mj/^ +1 ) dxidx 2 di ' 

da r , ( x\x 2 i \ d 3 



f dx d I X I (2 32) 

dP T J Xl V \x iy ^P T J d Xl dx 2 dt ' 



7 



And the allowed regions of the variables are given by 

s 

s = M 2 j^ +1 = x x x 2 s > M 2 ^, 

0<x T < ^ 2 ~ r \ (2.33) 

y/XiX 2 

2.3 The Nonrelativistic-QCD (NRQCD) Factorization For- 
malism 

First we consider the general method to get the NRQCD cross section for the process 
a + b — > (QQ)( 2S+1 Lj )(— > H) + c, where H is the final state heavy quarkonium and 
(QQ)( 2S+1 L ( }' 8) ) is the intermediate (QQ) pair which has the corresponding spectro- 
scopic state. From now on, we use the subscript n to represent the spectroscopic 
(QQ) state of ( 2S+1 Lj n_1 ' 8 ^), for simplicity. Once the on-shell scattering amplitude 
of the process A(a + b — > Q + Q + c) is given, we can expand the amplitude in 
terms of relative momentum q of the quarks inside the bound state because the 
quarks, which make up the bound state, are heavy. For more details of the method 
to deal with the heavy quarkonium production following the BBL formalism [07 



we refer to Refs. p8[| . The heavy quarkonium H production cross section 

a(pi) + b(p 2 ) -> (QQ) n (P) + c(p 3 ) -> H + c + X is given by 

^ (a( Pl ) + 6(p 2 ) - (QQ) n (P) + c(p 3 ) - + c + X) 



44^1?. (2-34) 



(1< 1 ^|^ / (a(p 1 )+6(p2)^(gg)n(P) + c(p 3 ))| 2 . (2.35) 



where 

i 

di IQixs 
Here, A4' is the amplitude of the process 

a( Pl ) + b(P 2 ) - (QQ) n (P) + c(p 3 ), (2.36) 

which can be obtained by integrating the free particle amplitude over the relative 
momentum of the quark inside the intermediate state (QQ) n (P), after projecting 
appropriate spectroscopic state Clebsch- Gordon coefficients. The parameter C n is 
defined by 

_ f 2N C (color - singlet), , 
° n ~ \ N 2 — 1 (color - octet). {Z ^ n 

And (0\O^\0) is the non-perturbative matrix element representing the transition 

(QQ)( 2S+1 L^ y )^H. (2.38) 

Finally, J denotes the angular momentum of the intermediate state (gg)( 25+1 Lj 1 ' 8 ' ) ), 
not of the physical state H. For the case of color-singlet intermediate state, which 
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has the same spectroscopic configuration with H, we can relate the matrix elements 
to the radial wave-function of the bound state as 



(orao> 



•17T 



(S — wave) 



C n x(2J+l) [ A|^(0)| 2 (P-wave). 



(2.39) 



For example, if we consider the ip production via ( 3 S[ ), ^S^), ( 3 S^), ( 3 Pq 8 ^), 
and ^P^ 8 "*) intermediate states, then the partonic subprocess cross sections 
are given by 

da 



dt 



(octet) 



8M C \ dt 



dt 



+ (0|(^ 3 if ) )|0>x£^( 3 P 



do-' 



3 D (8h 



df 



da 

-j (singlet) 
dt 



l^(0)| 2 ^,3c(l) ); 



M c 4vr dt 
after imposing the heavy quark spin symmetry 

(O|0 J /^( 3 pj 8 >)|O) = (2J + l)(0|O J/ ^( 3 P (8) )|0). 



(2.40) 



(2A1) 



The subprocess cross section for the color-singlet gluon-gluon fusion |13| is well 
known: 



da . A/i 
—j (singlet) = -— - 



s 2 (s - 4M C 2 ) 2 + P(t - 4M 2 ) 2 + m 2 (m - 4M ( 



(s - 4M 2 ) 2 (t - 4M 2 ) 2 (w - 4M 2 ) 2 
where the overall normalization Hi is defined as 

J\fi = - (4na s ) 2 (47ia)e 2 c M c 3 Gi(J/V>)- 
9 

The parameter Gi(J/ip), which is defined in NRQCD as 



(2.42) 



(2.43) 



3M 2 



Gi(J/V) 



2vrM 2 



l^5(0)| 2 , 



(2.44) 



is proportional to the probability of a color-singlet (cc) pair in the ( 3 S'j 1 ' ) ) partial wave 
state to form a physical J/ip state. It is related to the leptonic decay width 



r(J/V> -> / + n) = - vre^a 2 Gi(J/^), 

o 



(2.45) 



where e c = 2/3. From the measured leptonic decay rate of J/ip, one can extract 

GxiJ/ip) « 106 MeV. (2.46) 
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After including the radiative corrections of 0(a s ) with a s (M c ) = 0.27, this value is 
increased to ~ 184 MeV. Relativistic corrections tend to increase G\{J /ip) further to 
~ 195 MeV IE7I. 



For the case of color-octet subprocesses, we can use the average-squared am- 
plitude of the processes 7 + g (or q) — ► (cc)( 2S+1 L^) + g (or q) in Ref. f28]] , after 
crossing (k — > —k and qi — > —52) 

ElA-lfO? + <? - (cc)( 2S+l Lf) + 7) (5, f , u) 

= \\M'\ 2 { 1 + g^ (cc)C s+1 Lf)+g)(i,s,u), (2.47) 
El-M'i 2 (9 + q - (cc)( 2S+1 Lf) + 7) (s, t, u) 

= \\M'\ 2 ( 1 + q^ (cc)C S+1 Lf) + q)(t,s,u). (2.48) 

As previously explained, there is only one color-singlet subprocess, but 9 color- 
octet subprocesses are contributing to the process 

P + P (p) -> J/^ + l + X. (2.49) 

Whereas the initial partons for the color-singlet process are only the gluons, quarks 
can also be the initial partons for the color-octet subprocesses. For the ( 3 S[ 8 ^ ) channel, 
the gluon contribution corresponding to Fig. 4 (a) is absent because the effective 
vertex corresponding to Fig. 2 is vanishing. Using the parton level differential cross 
sections for various channels, in the next section we discuss in detail the p T and 
Mj/^ +1 distributions in the hadronic J/ip + 7 production. We also compare the 
color-octet contributions with the color-singlet gluon-gluon fusion contribution. 



2.4 NRQCD Results 

[For all the numerical results including the comparison of the color-singlet contribu- 
tion with the color-octet contributions for J /if) + 7 productions at the fixed target 
energy and Tevatron energy, see Ref. [fBJ]]. To conclude of this Section, we have 
considered the color-octet contributions to the associate J /if) + 7 production in the 
hadronic collisions, also compared to the conventional color-singlet gluon-gluon fusion 
contribution. Within the color-singlet model the J /if) + 7 production in the leading 
order is possible only through gluon-gluon fusion process. As we expected according 
to the naive power counting, we found that the color-octet contributions are signifi- 
cantly suppressed compared to the color-singlet gluon-gluon fusion process. For the 
case of the hadroproduction at the Tevatron energies, we found that there is a cross- 
over point at high p T where the color-octet contributions become dominant, and it 
is still larger than the fragmentation contributions shown in Ref. ]4^| . This sup- 
pression of the color-octet contributions in the associated production is qualitatively 
consistent with the recent result of Cacciari and Kramer J34J. They investigated the 
J /if) + 7 production via resolved photon collisions in HERA, and found the strong 
suppression on the color-octet contributions compared to the color-singlet one as at 
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3 7P 



100 GeV for p T > 1 GeV p 



res®) 



15% 



o . 



(2.50) 



Though most events at the Tevatron collider are in the region where the color-singlet 
gluon-gluon fusion contribution dominates, one could require an additional cut, for 
example: p T < 6 GeV, to guarantee that the color-singlet gluon-gluon fusion process 
remains one of the cleanest channels to probe the unpolarized [[C2| and polarized gluon 



distribution inside proton and to study heavy quarkonia production mechanism 



41 



In next two sections, strictly within the color-singlet model we consider the 
productions of J/ip + 7 at e + p and p + p(p) colliders to probe the unpolarized gluon 
densities (Section 3), as well as at the polarized fixed target experiment and polarized 
p + p collider to probe the polarized gluon densities of proton (Section 4). 



3 Probibg Unpolarized Gluon Densities from e + p 
and p + p(p) Colliders 

In this Section, we use the color-singlet model and leading order expressions every- 
where; we will also stick to leading order structure functions when evaluating ( |1.5| ). 
The NMC collaboration found || that this leading order formalism describes the 
shape of all kinematical distributions quite well, while the prediction for the normal- 
ization of the signal was too small by a factor of 2.4, even though the QCD corrected 
version of ( |1.7| ) was used, which increases the cross section by roughly 45%. This large 
"k-factor" is probably mostly due to the nonrelativistic treatment of the J/ip, in which 
all information about the wave function is contained in |i?(0)| 2 . One can therefore 
expect a similar k-factor for our reaction, so that our results for total cross sections 
should be considered as conservative estimates. In accord with our philosophy, we use 
the one-loop expression for a s with Nf = 4 active flavors and Aq C d = 200 MeV, and 
we take Q 2 = m 2 + p 2 T as momentum scale both in a s and in the structure functions. 

In a hadronic environment, a J/ip can probably only be identified when it decays 
into an e + e~ or pair; all our results therefore include a factor of 0.14, which is 

pSfl the combined branching ratio for these decays. Our final state thus consists of 
a hard, isolated photon in one hemisphere, balanced by an e + e~ or yU + /i~ pair in the 
opposite hemisphere, without any hadronic activity (other than the usual spectator 
jets resulting from the break-up of the incoming hadrons). This signal should be 
virtually free of any physical or instrumental backgrounds. We are now in a position 
to present predictions for J/tj) + 7 production at pp and ep colliders. 

We begin with a discussion of J/ip + 7 production in pp collisions. We focus 
on the Tevatron, which offers the largest cross sections of all existing colliders, and is 
expected to eventually accumulate an integrated luminosity of at least several hundred 
pb^ 1 . As discussed above, we require the J/ip meson to decay into a pair of charged 
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leptons. We then apply the following cuts, which should guarantee that the events 
are contained in the detectors and can be triggered upon: 



VI=V% >$GeV; (3.51) 
\y^\ < 3.5. (3.52) 

At present, the CDF detector can only detect muons with < 0.7, but future 
upgrades, as well as the upcoming DO detector, should provide better coverage; ( |3.52|) 
roughly describes the coverage of the electromagnetic calorimeters at CDF. 

In Figs. 6a,b we present the transverse momentum and energy spectra of the 
photon and the two leptons after the cuts (|3.51| , [3.52} ) have been applied, where we 



have used EHLQ1 J49| structure functions.^ In both figures we show the spectrum of 
the "harder" (denoted by 'b') and "softer" (denoted by 's') lepton separately, where 
the "hardness" is defined by the quantity plotted. Notice that the cut ( |3.51| ) implies 
that at least one lepton satisfies p l £ > 2.5 GeV, E lb > 2.9 GeV; together with the 
hard photon, this harder lepton can therefore be used to construct a trigger for this 
reaction. The transverse momentum and energy of the other lepton can in principle 
be arbitrarily small; however, Figs. 6 show that the additional cuts p l , E l > 
would not reduce the signal very much. On the other hand, it is necessary to include 
events where at least one lepton has less than 5 GeV transverse momentum. Due 
to the relatively mild cut (|3.52|) on the rapidities of the final state particles, the 
energy distributions are substantially harder than the p T spectra. Events where the 
whole J/tj) + 7 system undergoes a strong boost are interesting since they can yield 
information about the gluon densities at very small x, down to a few times 10~ 4 . 

In Fig. 7a we show the rapidity distribution of the produced photon. As 
discussed in Sec. 1, the overall normalization of the cross section is quite uncertain. 
We therefore normalize these distributions by dividing by the total cross section 
after cuts. In order to demonstrate the sensitivity of this distribution to the shape 
of the gluon density function f g \ p , we show results for the EHLQ1 ||9|] and D02 
|50| parametrizations. These two parametrizations make quite different assumptions 
about the large x behaviour of f g \ p : 

f g \p(x,Ql = 5 GeV 2 ) ocx- 1 (l + 3.5a;)(l -x) 5 - 9 (EHLQ1); (3.53) 
f g \ p (x, Ql = 4 GeV 2 ) ocx~ 1 (l + 9x)(l-x) 4 (D02). (3.54) 

The harder gluon distribution function of the D02 parametrization leads to a signif- 
icantly broader rapidity distribution; when going from y 1 = to |y 7 | = 3 the cross 
section only falls by a factor of 2.2, while the EHLQ1 gluon predicts a reduction by 
a factor of 2.9. 

The results of Fig. 7a have been obtained by integrating over p 7 and the rapid- 
ity of the J/ip. It is conceivable that two parametrizations of f g \ p which differ both 
in the large x and small x regions lead to similar results for the single differential 
cross section shown in this figure, since large photon rapidities correspond to very 



4 Using the recently available structure functions from MRS and CTEQ will not change the essence 
of physics given here, although it will certainly give different numerical results. 
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asymmetric initial states. Once a sufficiently large number of events has been ac- 
cumulated, such ambiguities can be resolved by studying the triple differential cross 
section da / (dpldy 1 dy^) . As an example, we show in Fig. 7b this triple differential 



cross section as a function of p^ 



at the symmetric point 



f 9lp (x ~ y/(2ml + 4(pi)*)/s) 
> 5 GeV then implies x > 6 



= 0. This cross 

where s is the 
10~ 3 , while the 
x > 0.02, even 



section is directly proportional to 
pp centre-of-mass energy; the cut p^ 

cross section becomes too small to be useful if p 1 > 15 GeV, i.e. 
if some 10 3 p6 _1 of data can be accumulated. As mentioned above, the range of x 
values that can be probed can be extended by studying more asymmetric configura- 
tions; it should therefore be quite easy to distinguish between parametrizations whose 
small x behaviour is governed by a different (negative) power of x, like the HMRS+/- 
parametrizations of Ref. | 51| , which assume x ■ f g \ p {x) oc x ±0 ' 5 . 

We now turn to a discussion of J/i/; + r y production at the upcoming ep collider 
HERA. Since in leading order this final state can only originate from a gg initial state, 
the observation of a sizeable signal would be an unambiguous proof for a nonvanishing 



gluon content of the photon. (The same is true [5^] for inclusive J / ip production in 77 
collisions.) In order to use (|1.5| ) for the cross section calculation, one has to convolute 
the gluon content of the photon with the photon content of the electron: 



where 



f"f\e(Z) Q ) 

We use the following cuts: 



Q 
7TZ 



(1 



hi 



Q 2 



p% > 1.5 GeV; 



-3.5 < yw < 3, 



(3.55) 



(3.56) 



(3.57) 
(3.58) 



where negative rapidities correspond to the proton beam direction. Notice that the 
cut ( |3.57[ ) is still sufficient to remove the contribution from reaction ( O ); the cut 
Q3.58j ) roughly describes the acceptance of the ZEUS detector. 

Notice that / fl | 7 is not constrained by a momentum sum rule; both the shape 
and the normalization of that function are unknown. It has recently been shown |53| 
that even a global fit to all existing data on F% does not yield much information on 
fg\ y . Our results should therefore be taken as examples, not QCD predictions. 

In Figs. 8a,b we show the transverse momentum and energy spectrum that 
result when the LAC1 parametrization |53| is used for f g \ T We see that the signal 



will only be useful if one lepton is allowed to have p T below 2 GeV; on the other 
hand the cut fl3.57| ) and the Jacobian peak in the p T distribution of the leptons 
relative to the axis of the J/ip imply that one lepton usually does have more than 
2 GeV transverse momentum. Also, due to the asymmetric nature of ep colliders, 
the difference between the energy distributions of Fig. 8a and the p T distributions 
of Fig. 8b is even more pronounced than for pp colliders; e.g., one could impose a 
cut E 1 > 5 GeV without loosing much signal. Such a cut might be helpful since the 
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resolution of electromagnetic calorimeters increases, and hence the relative errors on 
E y and p 1 shrink, with \[EP . 

In Fig. 9 we again compare the shape of the photon rapidity distributions as 
predicted by different parametrizations of photon structure functions. The present 
lack of data constraining / ff | 7 is reflected by the large differences between the three 
curves, corresponding to sets 1 and 3 of Ref. |53|] and the older DG parametrization 



of Ref. ||54[| , even though we have used a logarithmic scale for the |/-axis. In Ref. 
54] it was assumed that all gluons inside photons originate from radiation off quarks; 



in contrast, the analysis of Ref. |53| includes a truely "intrinsic" gluon inside the 
photon. At the input scale, the three parametrizations for / 9 | 7 have the following 
functional forms: 

f gh {x,Ql = 4GeV 2 ) ocx- L34 (l-x) 12 - 6 (LAC1); (3.59) 
f gh (x, Ql = 1 GeV 2 ) oc x 5 ' 9 (l - x) ' 56 (LAC3); (3.60) 
f gh {x, Ql = 1 GeV 2 ) oc x" L41 (l - a;) 4 - 5 (DG). (3.61) 

Clearly, the LAC3 parametrization is rather "pathological"; however, even this ex- 
tremely hard gluon distribution cannot be ruled out by data on alone. Fig. 9 
shows that it would lead to a rapidity distribution that is much more symmetric 
around y 1 = than the predictions of the other two parametrizations. Notice that 
the LAC2 gluon distribution falls even more rapidly at large x than is assumed in the 
LAC1 parametrization; it thus predicts even fewer events with y 1 > 0. We therefore 
conclude that a few dozen well measured J/ip + 7 events at HERA would suffice 
to distinguish between the DG and the three LAC parametrizations, using only the 
shape of the rapidity distribution. Of course, there is no guarantee that any of these 
parametrizations will describe the data. 



4 Probing Polarized Gluon Densities of Proton 

Interest in high energy spin physics has been recently revived with the result from 
(and interpretations thereof) the EMC collaboration []55J on polarized fi—p scattering. 
Processes in polarized pp collisions (such as achievable at an upgraded Fermilab fixed 
target facility or at a polarized collider |]56|] ) sensitive to the polarized gluon content 
of the proton, such as jets |57|, |59|, direct photons |60], |6"H , and heavy quark 



production f62[l , have been discussed. Another intriguing suggestion, due to Cortes 
and Pire is to consider X2{cc) production where the dominant lowest-order sub- 
process would be gg — > \2- The partonic level asymmetries for X2/X0 production have 
been calculated in the context of potential models p5| and are large. Low transverse 



momentum quarkonium production in polarized pp collisions using other methods has 
also been considered [^, [)6j as has high p T ip production p7fl . 

In all cases of charmonium production, the experimental signal is £ + £~ [I = e 
or fi) with the lepton-lepton invariant mass giving the J/ip mass, since xj can decay 
radiatively to J/ip + 7, and the J/ip signature is quite clean. As has been noted |j8] 



the question of extracting the gluon distribution is made less clean by the multitude 
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of contributing processes, e.g. 



9 + 9 ~ 


-> Xo,2 


9 + 9 ~ 


+ XJ + 9 


q + g - 


-> Xo,2 + q 


q + q - 


-* Xo,2 + 9 


9 + 9 ~ 


-> J/ip + g 


9 + 9 ~ 


-> &(— J/V> + X) + & 


q + q - 


-> 6(->J/V> + X)+&. 



(4.62) 



The simplicity of the Cortes and Pire idea is now gone. A full 0(a^) calculation of 
the spin-dependent production of xj is necessary. At low p T , \j production will also 
involve q + g and q + q initial states, while at high p T in addition the 2 — > 2 kinematics 
make the extraction of parton distribution functions less direct. Furthermore, a very 
careful calculation is required because even processes with small cross section can 
have a large effect on the asymmetry. The extraction of Ag(x,Q 2 ) using inclusive 
J /if) will be a challenge. 

Recently, J ftp produced in association with a 7 has been proposed as a clean 
channel to study the gluon distribution at hadron colliders ||12|| . The radiative xj 
decays can produce J/ip at both low and high p T , but the photon produced will be 
soft (E ~ (9(400 MeV)). If we insist that the experimental signature consist of a J ftp 
and 7, with large but equal and opposite p T there is only one production mechanism 
within the color-singlet model ||12||: 



g + g^Jji) + 1 . (4.63) 



Following Ref. p2| , this mechanism has been proposed in Ref. |)9j to study the po- 
larized gluon distribution in polarized fixed target experiments; we perform a more 
detailed analysis, including the analysis of this mechanism at the Brookhaven Rel- 
ativistic Heavy Ion Collider (RHIC) at both 50 GeV and 500 GeV center of mass 
energy and at the Large Hadron Collider (LHC). Polarized proton-proton operation 
is being considered for RHIC, for at least several months data collection, while the 
tunnel design of the LHC has been modified for the possible future inclusion of the 
Siberian Snakes needed for polarized proton-proton mode. Also, we list the full set 
of helicity amplitudes for this process, explicitly stating the Lorentz frame in which 
the J/ip helicities are given. 

The full helicity amplitudes for g + g — > J/ijj + 7 can be calculated following 
the approach of Gastmans, Troost and Wu [fTOH , with the addition of explicit helicity 
polarization vectors for the J/ip. A convenient set of polarization vectors can be 



found in Bohm and Sack ||71|| . These polarization vectors reduce to the usual massive 



vector boson (+,—,0) polarization vectors in the parton center of mass frame, and 
so, although the expressions for the helicity amplitudes have Lorentz invariant form, 
the (+,—,0) only refer to the J/ip helicity in this one particular frame. We find only 
one independent helicity amplitude (M(++, ++), where the '++,++' refer to the 
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helicity of gig2, iJ/ip respectively), and the remaining 5 non-zero helicity amplitudes 
can be found by crossing and parity symmetries: 

M(++,++) = M( , )=C 



s(s - 


- M 2 ) 




(s-M 2 )(t- 


M 2 )(u 


— M 2 ) 


u(u - 


- M 2 ) 




(s-M 2 )(t- 


M 2 )(u 


— M 2 ) 


i(i - 


- M 2 ) 




(s-M 2 )(t- 


M 2 ){u 


— M 2 ) 



M(+-,-+) = M(-+, +-) = C- — (1.0 1) 

Af (-+,-+) = M(+-,+-)=C- 



where C = ^ e <} e 9sR(,fyM5 — g ere m is the J lib mass, s, t and u are the usual 

Mandelstam variables, R(0) is the radial wavefunction at the origin of the cc in the 
.1 1 ib and a, b are the color indices of the incident gluons. Thus, the (spin and color) 
summed and averaged matrix element squared can be found 1]| : 



\M(g + g^J/ib + 1 )\ 2 = W«*M\m 



(t-M 2 ) 2 (u-M 2 ) 2 



^ (u-M 2 ) 2 (s-M 2 ) 2 ^ (s-M 2 ) 2 (t-M 2 ) 2 



(4.65) 



|i?(0)| can be related to the leptonic width of the J/ib: 

r(J/^e+e-) = ^|i?(0)| 2 = 4.72keV 

|i?(0)| 2 = 0.48 GeV 3 . (4.66) 

We are interested in the longitudinal spin-spin asymmetry, defined as: 

a( ++ )-a( + -) 

a ll = . > — 7- — r 4.67 
a(++) + a(+-) 

where <r(++) (<r(H — )) is the cross section for the collision of 2 protons with the same 
(opposite) helicities. This can be calculated in the parton model, 

All& = J dxi dx 2 d LL a Ag(xi,Q 2 ) Ag(x 2 ,Q 2 ) (4.68) 



where a is the parton level cross section (related to \M\ 2 given earlier), Ag(x,Q 2 ) 
is the polarized gluon distribution in the proton (= (g + (x, Q 2 ) — g~(x,Q 2 )) where 
g + (x,Q 2 ) (g~(x,Q 2 )) is the distribution for gluons with the same (opposite) helicity 
as that of the proton) and an is the parton level asymmetry 

&, *(++)-*(+-) (4 69) 

Given the known helicity amplitudes for this process, the parton level asymmetry is 
simply 

s 2 (s - M 2 ) 2 - t 2 (t - M 2 ) 2 - u 2 (u - M 2 ) 2 
aLL ~ s 2 (s - M 2 ) 2 + t 2 (t - M 2 ) 2 + u 2 (u - M 2 ) 2 ' ^ ' ^ 
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Measurable quantities of interest are the p T distribution and the joint p T — y\— 
y 2 distribution with 2/1 = 2/2 — 0, where yi( 2 ) is the rapidity of the 7 (J /if)). In the 
latter case, both partons have the same Bjorken-x (which is a function of p T only). 
The corresponding asymmetries are given by: 



A 1 

LL 



a{++) - a{+- 
a(++) + a(+- 
da(++) _ da(- 



" da(++) da(+-) 
dp T ~<Ip~ T 

dcr(++) I rfa(+-) I 

dp T dy 1 dy 2 lyi=m=0 dp T dy l dy 2 ^ 1= y 2=0 



A 3 

LL 



da{++) , rfff(+-) I 

dp T dy x dy 2 lw=w=° + dp T dyidy 2 ^ yi=y2=0 



Note that A| L is proportional to [Ag(x(p T ), Q 2 )] 2 . Another interesting theoretical 
concept (though not measurable experimentally) is the average or 'resolving 
power'. It is defined in the following way 

(d LL )<J = J dx! dx 2 d LL a g(x 1 , Q 2 ) g(x 2 , Q 2 )- (4.72) 

As we wish to determine if a given experimental scenario can shed light on the size of 
the polarized gluon in the proton, we need, in addition to calculating the asymmetry, 
to estimate the experimental uncertainty in the asymmetry. We will approximate the 
uncertainty by the statistical uncertainty, since ratios of cross sections should be rela- 
tively free of systematic uncertainties. The statistical uncertainty in the measurement 
of an asymmetry is given by 8 A, where 



= (473) 



and N is the number of events. 

We examine this process in several different experimental settings. First, we 
consider an hypothetical fixed target experiment and to be specific, take the proton 
beam energy to be 800 GeV (such as would exist at the upgraded Fermilab fixed 
target facility). In order to estimate the luminosity possible at such an experiment, 
we must make some assumptions. First, the Main Injector at Fermilab can provide 
~ 10 14 (unpolarized protons)/sec, with a 65% duty cycle [J72J. We'll assume a one 
month run, at a much reduced proton rate (say, a factor of 100), combined with a small 
polarized gas {H 2 ) jet target (approximately 1 cm long). This will give, we think, a 
very conservative estimate of / Cdt = 50 pb~ . We place no cuts on the rapidity of 
the photon or J/if), nor on the p T of the photon or leptons. We find a cross section of 
approximately 200 pb, most of which is at low p T . The resolving power (or average 
d^i) is found to be about 28%. We use the polarized distributions of Bourrely, Guillet 
and Chiappetta W\ . They provide 2 sets of distributions, one with a large polarized 
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gluon distribution and small polarized strange quark distribution (we'll refer to it 
as the set BGCO) and one with a moderately large polarized gluon and moderately 
large polarized strange quark distribution (we'll refer to this set as BGC1). The p T 
distribution is shown in Figure 10a (in cross section) and in Figure 10b (in A 2 LL ). 
We were also interested the asymmetry A 3 LL , (technically, instead of taking y± and 
y2 derivatives, we bin the events in the usual way, displaying the contents of the bin 
with —0.1 < yi,y2 < 0.1). The results are shown in Figure 12a (distribution in cross 
section) and 3b (A LL vs. p T ). We present in Table 1 the total number of events 
expected (at all p T and yi^ consistent with our cuts) as well as the 'resolving power' 
and asymmetry A l LL and an estimate of the statistical uncertainty, SA LL . We also 
list the number of events in a single p T bin (p T given in the table caption), and A\ L 
and 8A 2 LL for that particular p T bin. Finally, we present the the number of events in 
the same p T bin, further restricting the events to lie within |y 1)2 | < 0.1, and the value 
of A\ L and &A\ L in the particular p T bin. These are representative results. Higher 
statistics can be obtained by the inclusion of all p T bins. 

At this point, we would like to further address the work of Ref. p5| . The large 
asymmetries shown are surprising, and in our opinion not correct. The parton level 
asymmetry, making the following replacements for t and u {i.e. working in the parton 
center of mass frame): 

i = -i(s-M 2 )(l-cos#) 

u = -i(s-M 2 )(l + cosfl) (4.74) 

reduces to 

1 - |[(1 + 6 cos 2 6 + cos 4 9) + 2^(1 - cos 4 9) + f^(l - cos 2 8) 2 ] 
aLL ~ 1 + |[(1 + 6 cos 2 9 + cos 4 6) + ^f(l - cos 4 0) + f^(l - cos 2 9) 2 ] ' U '" )) 

Here cos 8 is measured in the parton center of mass frame. It is obvious that for 
cos9 = ±1, clll is a minimum (actually zero), and so, for any s, the maximum of cill 
should be at cos# = 0. In this limit, the asymmetry reduces to 

1 i( s + M 2 \ 2 

a LL (cos6 = 0) = 8 ^ 8 J 2 . (4.76) 

1 .1 f s + M 2 ^ 

Two further limiting cases are possible, namely production at threshold (s = M 2 ) 
which gives olll — | an d production at very high energy (s — > oo) which gives 
o-ll = §• For = — 38.75 GeV (the fixed target energy considered both here 
and in Ref. [§^]), the parton level asymmetry is near it's maximum value. Since 
Ag(x,Q 2 )/g(x,Q 2 ) < 1 generally, the maximum observable asymmetry is bounded 
by the maximum parton level asymmetry. Thus we are unable to understand the 
prediction, in Ref. |p5| , that the observable asymmetry can be as large as 85%. 

Next, we consider collider experiments at RHIC. RHIC is a high luminosity 
(£ = 2 x 10 32 cm~ 2 sec _1 = 6000 pb _1 /yr) collider capable of producing proton on 



18 



proton collisions for center of mass energies between 50 and 500 GeV. A program of 
polarized proton on proton collisions, at full energy and luminosity, is being discussed 
|T^| . We will assume a nominal running time of 2 months, at full luminosity, for 
50 GeV and 500 GeV each. In order to be somewhat conservative, we will estimate 
event numbers based on 300 pb _1 integrated luminosity. We will assume a generic 
collider type detector, and in order to simulate the acceptance we will require the 
photon and electrons observed to lie in the rapidity range \y\ < 2 (this simulates 
the acceptance of the proposed STAR detector at RHIC |75]| , level 2 for photons and 
electrons. We will not consider the possibility of the detection of the fi + final state 
at RHIC). Furthermore, we will (rather arbitrarily) require the p T of the photon larger 
than 1 GeV in the following discussion. We present our results for the p T distribution 
in Figure 12a, and A\ L in Figures 12b {yfi = 50 GeV) and 12c (^s = 500 GeV). 
See Figure 13a for ^ dy dy vs - Pt an d Figures 13b (-y/i = 50 GeV) and 13c 

(t/s = 500 GeV) for A\ L vs. p T . The 'resolving power' increases with energy (actually 
p T ), even though the observed asymmetry decreases. This is simply a consequence 
of the behavior of the polarized gluon distribution. Please refer to Table 1 for some 
representative results. 

In conclusion of this Section, we have studied the process p + p — > J Jijj + 7 + X 
in polarized proton-proton collisions. We first presented the necessary helicity ampli- 
tudes and discussed the calculation. Then we studied this process at polarized fixed 
target and in colliders, at polarized RHIC (50 and 500 GeV center of mass energy). 
Our results indicate that a polarized (double spin) fixed target program can be very 
useful in the determination of the polarized gluon distribution. It is unfortunate 
that no such experiment is planned. RHIC (especially at lower energies) is an ex- 
cellent probe of the polarized gluon distribution. Since A\ L is directly proportional 
to [Ag(x(p T ) , Q 2 ) / g(x(p T ) , Q 2 )} 2 , this distribution provides an easy determination of 
the polarized gluon distribution at various x values. It will prove especially useful to 
measure this distribution at several center of mass energies. Even a measurement of 
A\ L can provide much useful information (though it is not clear whether the higher 
statistics involved in this measurement will outweight the cleanliness of the extraction 
of the polarized gluon distribution in a measurement of A 3 LL ). The LHC can probe 
a much lower x in this process, and since Ag(x,Q 2 )/ g(x,Q 2 ) <C 1 there is no mea- 
surable asymmetry. However, the 'resolving power' at LHC can be still very large, 
so the smallness of the asymmetry is purely a consequence of the small-x behavior of 
Ag(x, Q 2 ). Polarized LHC can still be a useful tool for the study of high energy spin 
properties of the proton by utilizing a subprocess that will probe larger x (e.g. heavy 
Higgs production). We should also point out that we have considered only the color 
singlet model of heavy quarkonium production in this paper. A similar analysis can 
be performed using local duality, if it is determined at HERA that this mechanism 
contributes to J/ip + 7 production P0 |. Some slight modifications will be required, 
namely the inclusion of charm in the proton (this effect should be small) and light qq 
fusion, and in addition the modification of the parton level asymmetries. 
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Ntot 


{clll) 


Al L (SAi L ) 


n Pt 


A 2 LL (SAl L ) 


|yi, 2 | < o.i 


Al L (SAl L ) 


Fixed 


10500 


28.4% 


12.5% (1%) 


5000 


16% (1.4%) 


200 


22% (6%) 


Target 






3.2% (1%) 




4% (1.4%) 




5% (6%) 


RHIC 


11430 


43.3% 


19.1% (1%) 


4500 


26% (1.5%) 


1080 


32% (3%) 


50 GeV 






4.6% (1%) 




8% (1.5%) 




8% (3%) 


RHIC 


86400 


44.7% 


.4% (0.3%) 


4500 


1.7% (1.5%) 


840 


1.8% (3%) 


500 GeV 






.05% (0.3%) 




.2% (1.5%) 




.3% (3%) 



Table 1: Summary of representative predictions for J/ip + 7 production in polarized 
proton-proton interactions. Ntot is the total number of events above some minimum 
p T (= GeV for fixed target, 1 GeV for RHIC). (o-ll) is the 'resolving power' as 
defined in the text (this is independent of the polarized parton distributions). A\ L 
and 5A l LL are defined in the text; the upper entry corresponds to the large Ag(x, Q 2 ) 
(set BGC0) and the lower entry corresponds to the moderately large Ag(x, Q 2 ) (set 
BGC1). N p is the number of events in the particular p T bin (0.5-1.5 GeV for fixed 
target, 1-2 GeV for RHIC at 50 GeV, 3-5 GeV for RHIC at 500 GeV. 
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Figure Captions 

Figure 1 - Feynman diagrams for the color-singlet subprocess for g + g — > 
(cc)( 3 S[ 1 ) + 7 and the color-octet subprocess for g + g — > (cc)( 1 S , q 8 ' ) or 3 -Pj 8 ' ) ) + 7. 

Figure 2 - Feynman diagrams for the effective 7 — g — (cc) ( 1 5 , q 8 \ or 3 -Pj 8 ' ) ) vertex. 

Figure 3 - Feynman diagram for the effective q — q — (cc)( 3 S'j 8 ^) vertex. 

Figure 4 - Feynman diagrams for the color-octet contributions to the subpro- 

cesses (a) g + g -> [cc^S^, or 3 Pj 8) )+ 7 and (b) q + q ^ (cc)C S^ 8) , or 3 Pj 8) ) + 7 . 

Figure 5 - Feynman diagrams for the color-octet contributions to the subprocess 
q + q^ (cc)( 3 S{ 8) ) +7. 

Figure 6 - The transverse momentum (a) and energy (b) spectrum of the photon 
and the two leptons from J /if) + 7 production at the Tevatron collider with subse- 
quent J /if) — > e + e~,fi + fi~ decay, after the cuts of ( |3.51| , |3.52"D have been applied. The 



subscripts "b" and "s" refer to the lepton with the bigger and smaller p T (in a) or 
energy (in b), respectively. 

Figure 7 - The normalized rapidity distribution (a) and transverse momentum 
spectrum at the point y 1 = y^ = (b) of the photon from J /if) + 7 production at the 
Tevatron collider. Results for different parmetrizations of the gluon content of the 



proton are compared: solid curves - set 1 of Ref. E{J; dashed curves - set 2 of Ref. [50 



Figure 8 - The transverse momentum (a) and energy (b) spectrum of the 
photon and the two leptons from J /if) + 7 production at HERA with subsequent 
J /if) — > e + e~,fi + fi~ decay, after the cuts of ( |3.57| , |3.58"D have been applied. The mean- 
ing of the subscripts "b" and "s" is as in Fig. 2. 

Figure 9 - The normalized rapidity distribution of the photon from J/ip + j pro- 
duction at HERA. The solid, short dashed and long dashed curves have been obtained 
using the parametrizations of Ref. and sets 1 and 3 of Ref. |k| , respectively, for 
the gluon content of the photon. 
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Figure 10 - p T distribution, vs. p T (a) and A\ L vs. p T (b) for large 

Ag(x, Q 2 ) (solid line) and for moderately large Ag(x, Q 2 ) (dashed line) at fixed tar- 
get. 

Figure 11 - ^ ^^J yi=j/2=0 vs. p T (a) and A\ L vs. p T (b) for large Ag(x, Q 2 ) 
(solid line) and moderately large Ag(x, Q 2 ) (dashed line) at fixed target. 

Figure 12 - p T distribution, -fj^- vs. p T (a) for RHIC at y/s = 500 GeV 

(solid line) and at ^fs = 50 GeV (dot-dashed line), and A 2 LL vs. p T for RHIC at 
y/s = 50 GeV (b) and at y/s = 500 GeV (c) for large Ag(x,Q 2 ) (solid line) and for 
moderately large Ag(x,Q 2 ) (dashed line). 

Figure 13 - ^ ^ ^ \ yi=y2=0 vs. p T (a) for RHIC at \/s = 500 GeV (solid line) 

and at y/s = 50 GeV (dot-dashed line), and A\ L vs. p T for RHIC at \/s = 50 GeV 
(b) and at y/s = 500 GeV (c) for large Ag(x, Q 2 ) (solid line) and moderately large 
Ag(x,Q 2 ) (dashed line). 



[*) For Figures 1-5, please look Figs. 1-5 of |hep-ph/ 96 10294 . 

'*) For Figures 6-9, please look Figs. 2 - 5 of Z. Phys. C53 (1992) 673. 

[*) For Figures 10 - 13, please look Figs. 1 - 4 of Phys. Rev. D49 (1994) 4463. 
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